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Abstract 

The canonical decomposition of a Lorentz algebra element into a 
sum of orthogonal simple (decomposable) Lorentz bivectors is defined 
and discussed. This decomposition on the Lie algebra level leads to 
a natural decomposition of a proper orthochronous Lorentz transfor- 
mation into a product of commuting Lorentz transformations, each 
of which is the exponential of a simple bivector. While this later re- 
sult is known, we present novel formulas that are independent of the 
form of the Lorentz metric chosen. As an application of our meth- 
ods, we obtain an alternative method of deriving the formulas for the 
exponential and logarithm for Lorentz transformations. 

1 Introduction 

It is common to define a three-dimensional rotation by specifying its axis 
and rotation angle. However, since the axis of rotation is fixed, we may also 
define the rotation in terms of the two-plane orthogonal to the axis. This 
point of view is useful in higher dimensions, where one no longer has a single 
rotation axis. Indeed we may define a rotation in n-dimensional space by 
choosing a two-plane and a rotation angle in that plane; vectors orthogonal 
to the plane are fixed by the rotation. The most general rotation, however, 
is a composition of such planar, or simple, rotations. Moreover, the planes 
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of rotation may be taken to be mutually orthogonal to each other, so that 
the simple factors commute with one another. 

On the Lie algebra level, where rotations correspond to antisymmetric 
matrices, a simple rotation corresponds to a decomposable matrix; i.e., a 
matrix of the form u A v = uv'^ — vu^, where u, v are (column) vectors. 
Geometrically, the vectors u, v span the plane of rotation. A general rotation 
corresponds to the sum of decomposable matrices: A = J^i'^i ^ 
the mutual orthogonality of the two-planes of rotation corresponds to the 
mutual annihilation of summands: {ui A Vi){uj A Vj) = for i 7^ j. Note 
that mutual annihilation trivially implies commutation, so that exp(^jMj A 
Vi) = Yli exp('Uj Af j), and we recover algebraically the commutation of simple 
factors in a rotation. 

In the particular case of four-dimensional Euclidean space, a general ro- 
tation is the composition of two simple commuting rotations in orthogonal 
two-planes. Correspondingly, a 4 x 4 antisymmetric matrix can be written 
as the sum of two mutually annihilating decomposable matrices. It is possi- 
ble (although apparently not of sufficient interest to be publishable) to write 
down explicit formulas for this decomposition — both on the Lie algebra level 
and on the level of rotations. It is natural to ask if such an orthogonal de- 
composition is possible for Lorentz transformations. This article answers this 
question in the affirmative. If one assumes a Minkowski metric and allows 
the use of complex numbers, a Wick rotation may be applied to the formulas 
for four-dimensional Euclidean rotations. However our discussion will not 
proceed along this line; instead, we will limit ourselves to algebra over the 
reals (see below) and develop the formulas from first principles, without any 
assumptions on the specific form of the Lorentz metric. 

It should be noted that the orthogonal decomposition of a Lorentz trans- 
formation A into commuting factors is distinct from the usual polar decom- 
position A = BR, where 5 is a boost and i? is a rotation. While B and R are 
simple in the sense that they are exponentials of decomposable Lie algebra 
elements: B = exp(6) and R = exp(r), they do not necessarily commute, 
nor do the Lie algebra elements b, r annihilate each other. In particular, 
exp(6 + r) ^ exp(6) exp(r) in general. 

From the point of view of the geometry of Lorentz transformations, there 
is nothing new presented in this work. It is known that any proper or- 
thochronous Lorentz transformation is the product of commuting factors 
which act only on orthogonal two-planes ([5]). That the corresponding 
Lorentz algebra element can be written as the sum of mutually annihilat- 
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ing simple elements is a natural consequence. However, the formulas given 
here, which actually allow one to compute these decompositions using only 
basic arithmetic operations, appear to be original. Moreover they are not 
metric-specific, and so are applicable to any space-time manifold in a natu- 
ral way, without the need to use Riemann normal coordinates. 

Since it may seem oddly unnecessary to restrict ourselves only to alge- 
bra over the reals, we briefly explain our motivation for this. The author 
is interested in efficient numerical interpolation of Lorentz transformations, 
with a view towards a real-time relativistic simulation. That is, we are given 
(proper) Lorentz transformations Aq and Ai, which we think of as giving co- 
ordinate transformations, measured at times t — and t — 1, between an ob- 
server and an object in nonuniform motion. For any time t with < t < 1, we 
would like to find a Lorentz transformation A(t) that is "between" Aq and Ai. 
Analogous to linear interpolation, we may use geodesic interpolation in the 
manifold of all Lorentz transformations: A{t) — Aq exp(t log Aq ^Ai). Thus 
we need to efficiently compute the logarithm and exponential for Lorentz 
transformations. Although many computer languages offer arithmetic over 
the complex numbers, not all do, and in any case they are typically not 
supported at the hardware level, so their usage introduces additional com- 
putational overhead. 

2 Preliminaries 
2.1 Generalized trace 

We will make use of the k—th order trace tr^M of a n x n matrix M, which 
is defined by the identity 



Equivalently, the {n — k)-th degree term of the characteristic polynomial 
det(A/ — M) of M has coefficient (— l)^trfcM. In particular, the zeroth order 
trace is unity, the first order trace is the usual trace, and the n-th order trace 
is the determinant. For k 1, tr^M is not linear. However, the identities 
(i) trfciif^ = trfcM, (n) tikMN = tikNM, and (ni) tr^aM = aHvkM for 
any scalar a, are always satisfied for any k. These properties follow from 
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equation ([T]) and the basic properties of determinants. Of primary interest 
is the second order trace, for which we have the exphcit formula 

tr2M= ^(tr^M-trM^) (2) 

This follows by computing the second derivative of equation ([T]) and evalu- 
ating at t = 0. In engineering, tr2M is often denoted by IIm- 



2.2 Antisymmetry and decomposability 

In preparation for the next section, let us derive a few facts that we will need 
concerning 4x4 antisymmetric matrices. Such a matrix can be written in 
the form 

/ I ^ \ 

= ( -X w) ^^"^'"^ = h/3 -y^ (3) 
^ \-y2 yi / 

for some x, y G M.^. Here Wy is chosen so that it has the defining property 
WyV = y X V. 

By computing the determinant of (j3]) directly, we obtain the following, 
which is a special case of the general fact that the determinant of an anti- 
symmetric matrix (in any dimension) is the square of its Pfafiian. 

Lemma 2.1. det = (x ■ y)^ □ 

We may identify the wedge ( exterior) product of four-vectors u, v with the 
matrix u Av = uv'^ — vu^. This matrix is evidently antisymmetric, and if we 
write u = {uq, u), t> = (t>o, v), we find that uAv = Axy, where x = uqv — vqU 
and y = V X u. In particular, lemma 12.11 implies that det u Av = 0. This is 
again a special case of a general fact: the determinant of a wedge is zero in 
dimensions greater than two, since we can always find a vector orthogonal to 
the vectors u,v. In dimension four, the converse is also true. 

Lemma 2.2. A^y is a wedge-product if and only if det = 0. 

Proof. For the sufficiency, first assume that y 7^ 0. Choose vectors u, v such 
that V X u = y. Since x ■ y = (lemma [2.11) . we may write x = uqv — vqU 
for some scalars Uq, Vq. If y = 0, then take u = (1, 0) and f = (0, x). □ 

In general, the vector space of two-forms in n dimensions can be identified 
with the space of n x n antisymmetric matrices. In dimensions n > 4, not 
every two-form is decomposable; i.e., the wedge of two one-forms. Lemma 
gives a simple test for decomposability in dimension n = 4. 



4 



3 Lorentz algebra element decomposition 



Let g he a Lorentz inner product on M^; i.e., g is symmetric, nondegenerate, 
and det g < 0. Recall that the Lorentz algebra so{g) is the collection of 
all linear transformations L on for which L^g + gL = 0; or equivalently, 
= —gLg~^. We will refer to an element of so{g) as a bivector. 

3.1 Algebra invariants and properties 
Lemma 3.1. Ag G so{g) if and only if A is antisymmetric. 

Proof. Since g is symmetric, {Ag)^g + g{Ag) = g{A^ + A)g, which is zero if 
and only if = -A. □ 

Thus any Lorentz bivector L can be written in the form L = Ag for 
some antisymmetric matrix A. Although this identification gives us a vector 
space isomorphism between so{g) and the space of all antisymmetric 4x4 
matrices, the identification is not compatible with their respective Lie algebra 
structures; i.e., does not give a Lie algebra isomorphism. 

Lemma 3.2. For any L G so{g), tiL = = ti^L. 

Proof tikL = trfcL^ = tTk{-gLg-') = tik{-L) = (-l)HrfcL. □ 
Lemma 3.3. det L < for all L G so{g). 

Proof Write L = Ag. Then det L = det A det fi- < 0, by lemma O □ 

By the comments in section 12. the characteristic equation for a Lorentz 
bivector L is 

A^ + (tr2L)A2 + detL = 0. (4) 
The squared-roots (that is, the roots of + (tT2L)x + det L = 0) are 



At± = I (-traL ± y tr|L - 4 det L) (5) 

Since detL < 0, ^± are real numbers. Note that they are solutions to the 
equations fi+ + fi- = — tra-L and = det L, and that /x+ > and /i_ < 0. 

Moreover, det L 7^ if and only if /i+ > and yU_ < 0. 

Lemma 3.4. + (traL)^^ + (det L)I = 0. 
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Proof. This follows from equation and the Cayley-Hamilton theorem. □ 

The determinant and second order trace of a Lorentz bivector may be 
computed from the traces of its powers. 

Lemma 3.5. trsL = -^trL^ and det L = ^(tr^L^ - 2trL^). 

Proof. Lemma 13.21 and equation ([2]) imply the first identity. For the second, 
lemma El implies trL"^ + tr2L trL^ + 4 det L = 0. □ 

3.2 Simple bivectors 

Given two four-vectors u, v, we may construct the simple Lorentz bivector 
u V, which is the 4x4 matrix 

u V = uv^ g — vu^g (6) 

In covariant coordinate notation, {u A^ v)^ = u^Va — v^Ua- Observe that 
{u A^ v) = {u A v)g, whence lemma [3?T] implies that u A^ v & so{g). 

Lemma 3.6. L G so{g) is simple if and only if det L = 0. 

Proof. Write L = Ag, with A antisymmetric. Since detL = detAdetg, the 
lemma follows from lemma 12. 2[ □ 

Lemma 3.7. L G so{g) is simple if and only if + (tr2L)L = 0. Moreover 
if L = u A^ V, then tr2-L = [v^ gu){y^ gv) — {vFgvY . 

Proof. For L = u A^ v, one computes directly from ([6]) that (*) trL^ = 27 
and = 7L, where 7 is the scalar {vFgvY — {u^ gu){v^ gv) . Lemma [3.51 and 
(*) imply that ti2L = —7. For the converse, + {tT2L)L = implies that 
L^ + (tr2L)L2 = 0. Lemma [33] then yields detL = -\ltrL^ + tT2LtrL'^) = 0, 
and lemma 13.61 applies. □ 

The four-vectors u, v are parallel if and only if m A^f = 0. However, even 
if u,v are linearly independent (so that m f 7^ 0), it is possible that the 
metric is no longer nondegenerate when restricted to the plane spanned by 
u,v. We will call such a two-plane degenerate. The next lemma (combined 
with the previous lemma) shows that degenerate two-planes are synonymous 
with null two-planes (two-flats) in the sense of [6J. 

Lemma 3.8. u, v span a nondegenerate two-plane if and only ifti2uA^v 7^ 0. 
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Proof. For w = au + f3v, we have w'^gu = a{u^gu) + f3{v'^gu) and w'^gv = 
a{u^ gv) + (3{v'^gv). This hnear system has a unique solution if and only if 
its determinant {u^ gu){v'^ gv) — [u^gvY = tr2L is nonzero. □ 

The following lemma implies that almost all simple Lorentz bivectors are 
determined up to constant multiple by their squares. 

Lemma 3.9. Suppose L = uA^v is such that tr2L 7^ 0. Then Pl = — L^/tr2l/ 
is orthogonal projection ( with respect to g) onto the nondegenerate two-plane 
spanned by u,v; that is, Pi = Pl, trP^ = 2, and {gPiY = qPl- Conversely, 
if P is orthogonal projection onto a two-plane, the two-plane is necessarily 
nondegenerate, and we have P = Pl, where L = u v and u,v are any 
linearly independent four-vectors in the image of P. 

Proof. The first statement follows from lemmas [3171 13. 5[ and l3.8[ For the sec- 
ond statement, let V denote the two-plane forming the image of P. Suppose 
that w & V is such that w^gy = for all y & V. Thus for any four-vector 
X, w^gx = {Pw)^gx = w'^gPx = 0; hence w = 0, since g is nondegenerate. 
Thus V is nondegenerate. The remainder of the lemma follows from the 
uniqueness of orthogonal projection. Indeed, given a four-vector x, for any 
y E V we have {Px)^gy = x^gPy = x^gy. Since g is nondegenerate on V, 
the value of Px is thus uniquely determined by x and g. □ 



3.3 Orthogonal sum of simple bivectors 

Although an element L G so{g) will not be a simple bivector in general, we 
may write it as a sum of simple bivectors. Such a sum is not unique; however, 
we will do so in a canonical way using the operators 

= (7) 

Here fi± are defined as in equation (|5]). For P± to be well-defined, L cannot 
be a simple bivector; that is, det L ^ 0, since by the comments following 
equation (E]), this guarantees that /i+ — /i_ 7^ 0. 

Theorem 1. If L E so{g) is such that det L 7^ 0, then (i) P^ + P_ = I, 
(ii) P+P_ = = P-P+, (ill) Pl = P±, and (iv) P±L = LP±. Moreover, 
P±L are simple Lorentz bivectors with L = P+L + P-L, {Pj^L){P^L) = = 
{P_L){P+L), and trgPiL = -/i±. 
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Proof. Set AT = /i+ - Then + (traL)^^ + (det L)I = {L^ - I2+I){L^ - 
= -N^P^P_. Thus P+P_ = by lemma [131 Properties (i), (iii), 
(iv) also follow from the definitions of P±, coupled with lemma and the 
remarks after equation (|5]). The remaining two algebraic statements follow 
from (i), (ii), and (iv), so it suffices to establish that P±L are simple Lorentz 
bivectors with the second order trace as stated. First, note that (L^)'^ = 
{L'^)^ = {—gLg~^Y = ~gL^9^^- It then follows from equation ([7]) that 
\p±LY = -g{P±L)g~^; hence P±L G so{g). Second, {P±L)^ = P±L^ = 
±{L^ - jj,^L^)/N. Using lemma |3]11 we have = -(tr2L)L^ - (det L)L. 
Therefore, {P±Lf = ±fi±{L^ - l^^L)/N = n±P±L. From lemma |321 P±L 
is a simple bivector with ii2P±L = —fi±. □ 

Corollary 3.1. Any nonsimple Lorentz bivector L may be written as the 
sum of two simple Lorentz bivectors: L = + L_ with L^L_ = = L_L^, 
where 

and tT2L± = —fi±, with fi± are defined as in equation (jS]). □ 

We shall refer to the decomposition L = L+ + L_ in the corollary as 
the orthogonal decomposition of L. We remark that since /i+ > and 
/i_ < for a nonsimple Lorentz bivector, lemma 13.71 implies that the two- 
plane associated to L+ (as in lemma 13. 9p is time-like (intersects the null- 
cone) while that of L_ is space-like (does not intersect the null-cone) in 
Synge's classification of two-planes given in [6]. 



3.4 Special case: Minkowski metric 

In the case when g = rj = diag(— 1, 1, 1, 1) is the Minkowski metric, we may 
make use of the explicit parametrization of antisymmetric matrices given in 
equation (|3]). That is, any element of 50(77) can be written in the form 

= (x W^) 

for some x, y G M"^; i.e., Lxy = ^xy'?- Note that det Lxy = — (x ■ y)^, so that 
Lxy is simple if and only if x ■ y = 0. One computes, using lemma [331 that 
traLxy = |yp - |x|2. 
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Theorem 2. // x ■ y 7^ 0, then Lxy = i^a+b+ + -^a_b_ is the orthogonal 
decomposition of L^y: i^a+b+-^a_b_ = = iva_b_-^a+b+; where 

/i±x+(x-y)y u±y-(x-y)x 
a± = ± ana b± = ± 

and /i± = i(|x|2 - |y|2 ± ^(|x|2 - |y |2)2 + 4(x ■ y)2) . 

Proof. Using equation ([H]), compute L^y = i^x'y') where x' = (|xp — |yp)x + 
(x ■ y)y, and y' = — (x • y)x + (|xp — |yp)y. Now use coroUary 13.11 □ 

One may show that a_|_-a_ = 0, b_(_-b_ = 0, and a± xb^p = 0. Necessarily, 
a± ■ b-t = 0, since La+b+ and -La_b_ are simple. 

3.4.1 Relation with the Hodge dual 

We may also interpret orthogonal decomposition in terms of Hodge duals. 
Write the Lorentz bivector L as the two-form L = \L^e^ A e^, where 
(a = 0, 1,2,3) is a basis of and = r]°'^ea is the dual basis (the sum- 
mation convention is used). Observe that we may do this as Lri~^ = Lrj 
is antisymmetric. The Hodge dual (with respect to rj) of is A(_y)x, so 
that *Lxy = L(_y)x. Moreover, one computes that (x • y)a± = /i±b=p and 
(x ■ y)b± = — /i±azp. It then follows that 

*J^a+h+ = ^a-h- and * La_h- = ^a+h+- 

x-y /i+ 

3.4.2 Relation with the spin group 

In addition, we describe the orthogonal decomposition in terms of the double 
covering group homomorphism \E' : S'LgC — )■ SO'^l't]). Since we are using the 
mathematicians choice of metric signature, we use p{u) = {^^u^ ui-iu2) 
our linear embedding of into the algebra of 2 x 2 complex matrices, along 
with the involution {c d) ^^^^ P(^)* = P(^) and det p(u) = 

u^rju. The map ^ sends the 2x2 complex matrix M with unit determinant 
to the Lorentz transformation u 1— )■ [M p{u) M*) . 

On the Lie algebra level, the map ip : SI2C — )■ 50(77) corresponding to 
\l/ sends the traceless complex matrix m = ( " -a ) ^he Lorentz bivector 
u I— 7- p~^{mp{u) + p{u)m*). One computes that ip{rn) = Lxy, where x = 
{Re{b + c)Jm{b - c),2Re{a)) and y = {lm{b + c),-Re{b - c),2/m(a)). 
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Thus X ■ y = 2Im{a? + he) and |yp — |xp = —4:Re{a? + he). In particular, 
z/'(m) is simple if and only if + 6c is real. While it is possible to work 
out V'~^(-^a+b+), the resulting expression is not particularly nice. However, 
one computes that ip{im) = *L^y, so that by the previous discussion, if 
ipirrij^) = La+b+, then ip{i9m+) = La_b_ for some real 9. 

To get back to the Lie group level, we only need to exponentiate: if 
ilj{m) = L, then \E'(exp(m)) = exp(L). We will return to this briefly at the 
end of section 15.41 



4 Exponential on so{g) 



Recall that the exponential of a square matrix M is defined as exp(M) = 
Yl'^=o h.^^- "^^^ series always converges, and the exponential of a Lorentz 
bivector lies in the identity component the Lie group of all Lorentz transfor- 
mations. 

Theorem 3. If L is a simple Lorentz biveetor, then exp(L) = I + diL + d2LP' , 
where 



sinhv^^ti^ , , 1 - cosh A/^^tr^ r 

di = ^ — ana d2 = if tr2-L < U, 

V-tr2L tr2L 

sin ytr^ 1 - cos ^/ti^ 

di = — 1^=^- ana d2 = ^J tT2L > 0, 

Vtr2L tr2L 

and if tT2L = 0, then di = 1 and d2 = \- 

Proof. Lemma [32] implies that L^'^'^^ = (— l)^(tr2L)^-L for all A; 7^ 0, which 
is used to sum the exponential series. □ 

This formula appears in |2j with — tr2L replaced by the symbol a; however 
in that reference, a is not identified with a matrix invariant, the coefficient of 
LP' is incorrect, and the limiting case when ii2L = is not handled explicitly. 

For nonsimple Lorentz bivectors, one may sum the exponential series 
using the identity L^ = fi'^L^ + fi'^L_ for all k > 0, which one verifies by 
induction in conjunction with lemma 13.41 The resulting formula, which is 
cubic in L appears in |2] and in [1] , although the formula in the latter takes a 
different form than that in the former; both references use somewhat different 
methods than ours, and in particular, make use of algebra over the complex 
numbers. 
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Alternatively, we may make use of corollary 13 .11 to obtain the exponential 
of a nonsimple Lorentz algebra element. Since L+,L_ trivially commute, 
we have exp(L) = exp(L+) exp(L_), and each factor can be computed from 
theorem [3] to obtain the following. 

Theorem 4. If L E so{g) is nonsimple and L = + L_ is the orthogonal 
decomposition of L as in corollary then 



Note that if we write L± in terms of L (according to corollary 13. ip . we 
obtain a polynomial of degree six in L. However, the degree can be reduced 
to three by lemma 13. 4( the resulting expression is necessarily the same as 
that obtained in [1] and j2]. 

5 Lorentz transformation factorization 

The Lorentz group 0{g) is set of all linear transformations A on that 
preserve the metric: A'^gA = g. Necessarily det A = ±1, and A is said to be 
proper if det A = 1. It is well-known that the group of all proper Lorentz 
transformations has two connected components; the Lorentz transformations 
within the identity component are said to be orthochronous. In the case 
when g = 7] is the Minkowski metric, the orthochronous transformations are 
characterized by the component Aq being positive. The set of all proper 
orthochronous Lorentz group elements will be denoted by SO'^{g). 

A note on proofs. In the remainder, all proofs will use the following 
abbreviations without warning: = tr^A, = tr^L, x = a/— I2 when 
I2 < 0, y = y/l2 when I2 > 0, c+ = coshx, c_ = cosy, s+ = sinhx/x, and 
S- = siny/y. Note that since A~^ = g~^A^g, we have tr^A"^ = Tk = tr^A. 
We remind the reader that tr2A = |(tr^A — trA^), from equation ([2]). 

5.1 Some Lorentz transformation relations 



Lemma 5.1. For all A G SO^{g), ti^A = trA. Consequently, A satisfies the 
relation A^ - (trA)A3 + (tr2A)A2 - (trA) A + / = 0. 



exp(L) = / + d'lL+ + d^L^ + d^L^ + d^Li 
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Proof. The characteristic polynomials of A and A ^ are the same, namely 
— TiA^ + r2A^ — T3X + 1. By the Cayley-Hamilton theorem, (*) A"' — tiA^ + 
T2A'^-T3A + I = and (**) A-^'-riA-^ + rzA-^-rgA-i + J = 0. Muhiplying 
(**) by A^ and comparing with (*), we get T3 = ti. □ 

Lemma 5.2. If A e SO+{g), then A'^ + A-'^ = AJ + ^^(A + A-^), where 
A2 = -traA, B2 = tiA, A3 = {2 - tr2A)tr A, ^3 = tr^A - trzA - 1, A^ = 
(2 - traA) tr^A + tr^A - 2, and B4 = (tr^A - 2 trsA) trA 

Proof. From lemma [F!Tl we have A^ — tiA + T2I — tiA^^ + A^^ = 0; so that 
A2 + A-2 = -rgj + ri(A + A"^). Moreover, A^ - riA^ + T2A - rj + A^^ = 
and A-3 - riA-2 + rsA-^ - tJ + A = 0; thus A^ + A'^ = ri(A2 + A-^) - 
{t2 + 1)(A + A-i) + 2ri/ = ri(2 - ts)/ + (n^a - - 1)(A + A"^). Similarly, 
we have A^ + A"^ = ti{A^ + A'^) - T2{A^ + A-^) + ri(A + A-^) - 21 = 
(n^s - T2A2 - 2)1 + (nEg - T2B2 + ri)(A + A-i). □ 

5.2 Simple Lorentz transformations 

We will say that a Lorentz transformation is simple if it is the exponential 
of a simple Lorentz bivector. Our first goal is to give an algebraic criterion 
for simplicity. We will make use of the fact that the orthochronous Lorentz 
group is exponential (see [3]); that is, exp : so{g) SO^{g) is surjective. 

Lemma 5.3. For nonsimple A G SO^{g), there is a nonsimple Lorentz 
bivector L with orthogonal decomposition L = + L_ such that A = exp(L) 
andii2L± = — /i±, where y/Ji+ = cosh~^ i(trA+v^), -\/— /U_ = cos^"*^ ^(trA— 
\/A), and A = tr^A — 4tr2A + 8. Moreover, we have trA = 2(c+ + c_) and 
tr2A = 4c+c_ + 2, where c+ = cosh ^/J^^, c_ = cos ^—ji^. 

Proof. Since A is not simple, A = exp(L) for some nonsimple Lorentz bivector 
L, and theorem H] applies. Moreover, using lemma [3771 one computes 

A^ = / + 2rf+(l + fx+dt)L+ + 2rf^(l + fx_d2)L^ + A+Ll + A_Ll (9) 

where A± = {2df + d^"^ + f^±df'^) = 2(c^ — l)//i±. Computing the traces 
in theorem m and equation (Q, we obtain trA = 4 + c/^trL^ + (ig^^-^- 
trA^ = 4 + A+trLl + A_trL^_. Using trL| = -2tr2L± = 2/i±, one then 
computes that (*) trA = 2(c+ + c_) and trA^ = 4(c^ + c^ ^1)- Consequently, 
{■^^) tr2A = |(tr^A — trA^) = 4c+c_ + 2. Equations (^) and (*T*r) can then be 
solved for c± to obtain the formulas in the statement of the lemma. □ 
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Lemma 5.4. A is simple if and only if ii 2^ = 2trA — 2 and trA > 0. 

Proof. If A is simple, then A = exp(L) for some simple Lorentz bivector 
L. From theorem [3] and lemma 13. 5[ (*) ri = 4 — 2(i2^2- Moreover, we 
have A^ = J + 2diL + (2^2 + (i^)-L^ + 2did2L^ + c?2-Z^^, and one computes 
T2 = ^{Tf — trA^) = 6 + (df — 6d2)l2 + d\ll using lemma [33] and det L = 0. 
Thus (**) ra - 2ri + 2 = (rf? - 2^2)^2 + dill- By theorem El if Z2 < 0, 
then (i2^2 = 1 — cosh x and d\l2 = — sinh^ x = 1 — cosh^ x. Thus (*) yields 
Ti = 2(1 + coshx) > 4, and (**) yields r2 - 2ri + 2 = 0. Similarly, if I2 > 0, 
then Ti = 2(1 + cosy) > and T2 — 2ti + 2 = 0; and if I2 = 0, then ri = 4 
and r2 - 2ri + 2 = 0. 

Now suppose that A is not simple, and write A = exp(L) as in lemma 
15.31 Thus T2 — 2ti + 2 = — 4(c+ — 1)(1 — c_), which is only zero when c+ = 1 
or c_ = 1. As A is nonsimple, /i+ > and /i_ < 0; whence c+ > 1, and 
c_ = 1 only if y/—fIZ is a nonzero multiple of 2tt. However in the latter case, 
theorem m implies that exp(L) = exp(L+), which cannot be the case, since 
A is assumed to be nonsimple. □ 

We note that the proof of lemma 15.41 implies that if A is nonsimple, then 
c+ > 1 and — 1 < c_ < 1. In particular, trA > 0. 



5.3 Simple logarithm 

Theorem 5. Suppose A G SO~^{g) is simple. For trA > 0, let us define 
La = \k{A - A^i), where 

(i) k = — and ./a = cosh~VkrA - 1), trA > 4 
^ ^ sinh^/I ^2 ^' J 

(a) k = and = cos"^(itrA - 1), z/0 < trA < 4 

sm V — /i 

(Hi) k = 1 and ^ = 0, z/trA = 4 

Then La is a simple Lorentz bivector with exp(LA) = A and tra^A = — /i. 
For trA = 0, Pa = ^{I — A) is orthogonal projection onto a nondegenerate 
two-plane such that —tt^Pa is the square of a simple Lorentz bivector L with 
expL = A and tr2L = n^. 
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Thus if trA 7^ 0, then La may be taken as the logarithm of A: La = log A. 
In the case trA = 0, we may reconstruct log A from the two independent 
four-vectors in the image of Pa; it should be noted that in this case, A is 
necessarily an involution: A^ = /. In general, the logarithm is not unique, 
see In particular if L is simple with tr2L = 2mT, for any integer n, then 
A = exp(L) = / by theorem [31 In this case, the above theorem gives La = 0. 
A version of the formula for log A in the case trA 7^ that involves algebra 
over the complex numbers appears in |T]. 

Proof. Write A = exp(L) for some simple L. By theorem [31 we have (^) 
diL = |(A — A~^). As in the proof of lemma [531 (^) ri = 4 — 2(i2^2- 

From theorem [31 we see that di = only if > and y% is a nonzero 
multiple of tt; by the periodicity of the sine and cosine, we may assume that 
/2 = vr^. In this case, d2 = IjvP' and by t\ = 0. Conversely, if ri = 0, 
then (irk) implies (^2/2 = 2, which only happens if /2 > and is a nonzero 
multiple of tt. Thus if n = 0, then A = / + {2/tt^)L^ = I - 2Pl, with Pl as 
in lemma [331 Thus, Pl = Pa- 

If Ti 7^ 0, then di 7^ 0; and from (*) and theorem [31 we only need to 
deduce the value of I2 in terms of A in order to obtain the formulas (i) — (iii). 
If we assume that I2 < 0, then theorem [31 states that ^2^2 = 1 ~ coshx; and 
so (ick) can be solved to yield coshx = |ri — 1. Note that this equation has 
a solution for x > if and only if ti > 4. The cases when I2 > and I2 = 
are handled similarly. □ 



5.4 Decomposition into simple factors 

In the case when A is nonsimple, then L = log A is a nonsimple bivector, 
and we have the orthogonal decomposition L = L_(_ + L_ into simple mu- 
tually annihilating summands. Thus A = exp(L_|_) exp(L_) is a product of 
commuting simple factors. We give explicit formulas. 

Theorem 6. If A E SO^{g) is nonsimple, then A = exp(L) where L is a 
nonsimple Lorentz bivector with projection operators 

c+ — c_ 

where c± are as in lemma lSTR 
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Proof. From theorem HI we have (*) = exp(— L) = J — dfL^ — diL_ + 
d^L^ + (i^L^, and we obtain the equation |(A + A^^) = / + d^L"]^ + c?2 
which we rewrite as (o) c/^L^ + ^2 — ~^ + + ^ "'^)- addition, 
from equation (|9]), we have |(A^ + A"^) = / + A+L\ + y4_L^. However, 
by lemma I5.2[ this equation may be rewritten as (oo) A+Ll + A^L^_ = 



— \{t2 + 2)/ + |ti(A + A^^). Equations (o) and (oo) form a hnear system in 
L^, whose determinant is computed to be 2(c+ — 1)(1 — c_)(c_|_ — c^)/ fi+fi_, 
which is never zero. After inverting the system (o), (oo) , one computes 

= Ll + Ll = (;u_c+ - yu+c_)/ + i(^+ - /i_)(A + A"^) 

c+ — c_ 

Equation ([7]) then yields the desired formulas for P±. □ 

Theorem 7. Suppose A G 5*0+ ((7) zs nonsimple. Let c±, and fi± be defined 
as in lemma lSTSi If H- 7^ — vr^, then A = exp(L) with L a nonsimple Lorentz 
bivector whose orthogonal decomposition L = + L_ is given by 

where s+ = sinh y/JI+/ y/JI+ and s_ = sin ^J^-fiZI ^—{i^. Moreover, tr2L± = 
-/i±. 

Proof. From equation (|9]) and equation (*) in the proof of the previous theo- 
rem, we obtain the equations |(A — A^^) = dfL+ + d^L_ and |(A^ — A~^) = 
df{l + /i_|_(i^)L_|_ + di{l + yU„(i^)L_. These define a linear system for L±. 
The determinant is — (c+ — c_)s+s_, which is zero only when ^— is a 
nonzero multiple of vr. Inverting the linear system yields the formulas in the 
statement of the theorem. □ 

Theorem 8. For nonsimple A G SO^{g), A = A+A_, where A± are the 
commuting simple Lorentz transformations given by 

A± = ±7^ V {(1 + 2c±)/ - A-i - (1 + 2c^)A + A^} 

Z[C^ — C^) 
where c± are as in lemma lST^ 

Proof. Let L be a nonsimple Lorentz bivector such that A = exp(L), and P± 
as in theorem [6l Observe that exp(P±L) = Pzp + P±A. Indeed, {P±L)^ = 
P^L"' = P±L"' for all n > 0, since P± is a projection that commutes with 
L. Thus, exp(P±L) = I + ^^^^±(P^Lr = / + P±(E„>o ^^'^) = ^ + 
P-j-(exp(L) — /). The formulas for A± then follow from theorem O □ 
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As we noted in section I3.4.2[ in the special case when g = t] is the 
Minkowski metric, we may write L± in terms of the Lie algebra homomor- 
phism t/' : s/gC — ?■ 50(77). Namely, ?/'(m+) = L+ and ?/'(z^m+) = L_ for 
some traceless 2x2 complex matrix m+ = ( 2 -a ) such that + be is real. 
It follows that M+ = exp(m+) and M_ = exp(z6'm+) give (commuting) 
matrices in SL2C with \E'(M±) = A±; i.e., M± gives the decomposition of 
\I/(M+M„) = A into simple factors. Observe that exp(m_|_) is readily com- 
puted, since = (a^ + bc)I. 

We may view the decomposition of a nonsimple Lorentz transformation 
into commuting simple factors as a generalization of Synge's 4^screw: a 
product of a boost and a rotation in orthogonal two-planes. While it is 
known that every proper nonsimple Lorentz transformation can be expressed 
as the product of two commuting simple Lorentz transformations (see [5]), the 
formulas presented here are ostensibly original, and in any case independent 
of the specific form of the Lorentz metric. 
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